Antiferromagnetic fluctuations, symmetry and shape of the gap function 
in the electron-doped superconductors: the functional renormalization-group analysis 
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The problem of the symmetry of the superconducting pairing and the form of the gap function in the 
electron-doped superconductors is reconsidered within the temperature-cutoff functional renormal- 
ization group approach combined with the Bethe-Salpeter equations. The momentum dependence 
of the order parameter for antiferromagnetic and superconducting instabilities in these compounds 
is analyzed. The gap function in the antiferromagnetic (particle-hole) channel has its maxima at the 
hot-spots, or at the diagonal of the Brilloin zone in their absence. The wavefunction in the singlet 
superconducting channel is non-monotonic in the vicinity of the (tt, 0) and (0, 7r) points, deviating 
therefore from the conventional d-wave form in striking similarity with recent experimental data. 
An instability in the triplet superconducting channel is much weaker than the singlet one and has 
an f-wave like form of the gap function. 
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It is by now well established that the superconducting 
order parameter in high-T c compounds is described by a 
nearly {cosk x — cos fc y )-momentum dependence: its ab- 
solute value is largest at the Fermi surface (FS) points 
close to (tt, 0) and (0, tt) and vanishes at the FS crossings 
on the Brillouin zone (BZ) diagonals. 

Accurate measurements of the gap function in the hole- 
doped cuprates found, however, a slight deviation from 
this d x 2_ J/ 2-wave momentum dependence [1], with a flat- 
ter angular dependence near the nodal points. On the 
other hand, recent experiments on the electron-doped 
cuprates revealed non-monotonicity of the gap function 
as a function of an angle around the FS [2,3]. In particu- 
lar, the observed gap function has its maxima away from 
the points of the Fermi surface which are closest to the 
(tt, 0) and (0, tt) points of the Brillouin zone. Since the 
symmetry and the details of the momentum dependence 
of the superconducting order parameter are closely re- 
lated to the structure of the effective pairing interaction 
between the electrons, the momentum dependence of the 
gap function contains important information about the 
pairing mechanism in high-Tc compounds. 

One of the proposed candidates for the pairing mech- 
anism in cuprates is the antiferromagnetic (AFM) spin 
fluctuations scenario [4-7] . While the theoretical studies 
of the hole-doped compounds found that antiferromag- 
netic fluctuations lead to a d-wave symmetry of the gap 
function [8,9], there is currently a theoretical controversy 
about the symmetry of the superconducting pairing in 
the electron-doped cuprates. As it was argued earlier 
[2,3], with increasing electron doping the hot-spots of the 
Fermi surface (i.e. the points connected by the antifer- 
romagnetic wave vector Q = (tt, 7r)) move towards each 
other (see Fig. la). Since the sign of the d-wave gap is 
opposite at the hot-spots, this leads to the suppression 
of the tendency towards the d-wave pairing. 



In this situation, different types of the pairing sym- 
metry were proposed. An s-wave pairing at large elec- 
tron dopings was proposed [10], as arising mainly due 
to the effect of disorder in the model with both short- 
and long-range Coulomb parts of the interaction. On the 
other hand, for short-range interactions in the absence of 
disorder, p-wave superconductivity was considered as a 
possible candidate in the presence of strong charge fluctu- 
ations [11]. With these theoretical proposals, it remains 
an important problem which pairing symmetry is favored 
at intermediate and large electron dopings and what is 
the theoretically expected ground state superconducting 
order parameter for the electron-doped superconductors. 

To study this problem, we apply in the present paper 
the functional renormalization-group (fRG) technique 
[12-16], which proved successful in describing the inter- 
play of antiferromagnetism and d-wave superconductiv- 
ity in the t-t' Hubbard model. Although a previous ap- 
plication of this technique to the microscopic analysis 
of electron-doped superconductors within the Hubbard 
model found d-wave pairing even at relatively high dop- 
ings [17], the results of this study will be reconsidered in 
the present paper in two points. 

First, the momentum-cutoff renormalization group ap- 
proach of Ref. [17] does not allow to search for the charge 
and spin instabilities in the forward scattering (zero- 
momentum transfer) channel. This approach, therefore, 
may miss the possibility of the triplet pairing, which is of- 
ten closely related to such instabilities. This drawback is 
overcome in the temperature-cutoff fRG approach [15], 
which uses temperature as a natural cutoff parameter 
and proved successful in describing both AFM and fer- 
romagnetic instabilities together with singlet- and triplet 
superconducting pairing [15,16]. 

Second, a recent extension of the temperature cutoff 
renormalization group approach - its combination with 
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where ti 
and : 



: t for nearest neighbor (nn) sites i and j 
-t' for next-nn sites (t, t' > 0) on a square 
lattice. The chemical potential /i is determined by the 
filling n > 1; for convenience we have shifted the chemical 
potential fi by At' . 

The method. We follow the many-patch temperature 
cutoff fRG for one-particle irreducible Green functions 
as proposed in Ref. [15]. This scheme accounts for ex- 
citations with momenta far from and close to the FS, 
which is necessary for the description of instabilities aris- 
ing from zero-momentum transfer ph scattering. The mo- 
mentum dependence of the effective interaction between 
electrons at a given temperature, Vr(ki, k 2 , k 3 , k^) is 



FIG. 1. (a) The shape of the Fermi-surface (solid line) in 
the electron-doped superconductors, the dashed line shows 
the umklapp surface k x + k y — ir, black circles mark 
hot-spots, arrows indicate the direction of movement of the 
Fermi-surface and hot-spots with increasing doping, (b) The 
actual non-interacting Fermi surfaces of the model (1) at 
t' = 0.3t, and fj, = 0.8* (5 = n - 1 = 0.09) and fi = lAt 
(S — 0.26). The dots mark the positions of the centers of 
12 patches, located in the upper-right quarter of the Brilloin 
zone; a is an angle corrdinate of the patch centers. 



the Bethe-Salpeter (BS) approach - was proposed as a 
very convenient tool of studing symmetry and the form 
of the wavefunctions of different order parameters [18]. 
Contrary to previous fRG approaches, this method is not 
based on the knowledge of the wavefunctions at high tem- 
peratures. For the hole-doped superconductors it pre- 
dicts correctly the flattening of the gap function near the 
nodal points [18] in agreement with the experimental ob- 
servation [1]. 

In the present paper we use the temperature cutoff 
functional renormalization group approach in combina- 
tion with the Bethe-Salpeter equations to extract eigen- 
functions and eigenvalues of the effective interaction in 
the particle-particle (pp) or the particle-hole (ph) chan- 
nel in the electron-doped superconductors. 

The model. As a convenient model which describes 
electron-doped superconductors, we consider the 2D t-t' 
Hubbard model = H - (n - At')N with 



parametrized by the position of incoming ki , k 2 and out- 
going electron momentum k 3 at the Fermi surface. The 
fourth momentum, k4 is determined by the momentum 
conservation law. Neglecting the frequency dependence 
of the vertices, which is expected to have minor relevance 
in the weak-coupling regime, the RG differential equation 
for the interaction vertex has the form [15] 



dVr 



(2) 



dT dT dT 

where o is a short notation for summations over interme 
diate momenta and spin, 



Lp h ,pp(k,k') = 



Mgk) - M±gkQ 
£k =F £k' 



(3) 



and /r(£) is the Fermi function. The upper sign in 
Eq.(3) is for T p h and the lower sign for L pp , respec- 
tively. Eq.(2) has to be solved with the initial condition 
Vr (ki,k2,k3,k4) = U: the initial temperature is chosen 
as large as To = 400i. 

We discretize the momentum space in N p = 48 patches 
using the same patching scheme as in Ref. [15]. This re- 
duces the integro-differential equations (2) to a set of 
5824 differential equations, which were solved numeri- 
cally. The position of the centers of the patches at the 
Fermi surface for two different filling is shown in Fig. lb. 

To perform an analysis of possible instabilities within 
the fRG+BS approach, we consider the solution of the 
Bethe-Salpeter equations [18,19] 



i _ A ph <A 



£r£ h (k,p)£ ph (p,p + QX = T 



ph 



^r p n p (k,p)L pp (p,- P )0pp = T 
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The 2-particle reducible vertices Fp hpp (k, k') can be di- 
rectly extracted from the fRG flow according to 

f V T (k,k',k' + Q), ph (AFM) 

r^ P p(k,k') = I vHk,-k,k')±v r (k,-k,-k'), (5) 

[ pp (singlet, triplet SC) 

The value A p h,p P = 1 corresponds to an ordering in- 
stability with the symmetry of the eigenfunction ' pp . 
Therefore, tracing the temperature dependence of eigen- 
values and -functions allows to identify both, the lead- 
ing instabilities and their concomitant order parame- 
ter structure. We stop the fRG flow at the temper- 
ature Tx = O.OOlt (the maximum interaction vertex 
Vmax = max{V / (ki, k 2 ; k 3 , kj)} at this temperature re- 
mains smaller than the bandwidth). We have verified 
that the results for the eigenfunctions are only weakly 
dependent on the choice of Tx- 

Results. Below we discuss the results of the numerical 
solution of the Bethe-Salpeter equations in the electron- 
doping regime. We choose t' = 0.3i which is close to typ- 
ical values considered for the electron- doped cuprates. 
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FIG. 2. Eigenvalues (a) and angular-dependence on the FS 
of the eigenfunctions <j> p of the Bethe-Salpeter equation at 
T = Tx (b,c,d) for t' = 0.3i, U = 3.5t, fi = 0M (5 = 0.09). 
Tx is the lowest temperature reached in the fRG flow. The 
dots are the value of the eigenfunctions at the centers of cor- 
responding patches, cf. Fig. 1 



To remain in the weak-coupling regime, where the consid- 
ered technique is applicable, we put U — 3.5<. Although 
this value is substantially smaller than that which is ex- 
pected for the cuprate materials, we show that it allows 
us to reproduce the main features which are observed ex- 
perimentally. We also do not expect qualitative change 
of the results even for higher values of U. 

We start in Fig. 2 with the results for the chemical po- 
tential /i = 0.8i, which corresponds to a hlling n = 1.09, 
i.e. 9% of the electron doping. At this filling the eigen- 
value corresponding to the AFM instability saturates at 
a value A < 1 at the smallest temperature Tx which one 
can reach in the fRG flow. Therefore, the AFM is not 
expected to be the leading ground-state instability (more 
generally, the AFM is found not to dominate at fj, > 0.7t, 
i.e. S > 0.05). The corresponding eigenfunction in the 
particle-hole channel (Fig. 2b) has its maxima at the 
hot-spots of the Fermi surface. 

The eigenvalue A pp in the singlet pairing channel is 
monotonically increasing with decreasing temperature, 
and it is biggest close to Tx- Therefore, the singlet pair- 
ing is expected to be the leading instability at T — > 0. 
The corresponding wave function is non-monotonic (Fig. 
2c) and has a shape which is strikingly similar to the re- 
cent experimental data [2,3]. The maximum of the wave- 
function lies between the 3-rd and 4-th patch of the Fermi 
surface, i.e. it is shifted from the position of the hot-spot 
towards the point of the FS closest to the (n, 0) . The 
eigenvalue corresponding to the triplet pairing instabil- 
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FIG. 3. Same as Fig. 2 for fi = lAt (S = 0.26). 



ity is much smaller than for the singlet one, although it 
increases with decreasing temperature. The correspond- 
ing wave function (Fig. 2d) has nodes at the diagonals, 
and it is maximal near hot spots having therefore f-wave 
like symmetry. 

Now we consider the case of stronger doping /i = 1.40t, 
i.e. n = 1.26 (26% of the electron doping), see Fig. 3. 
In this case the hot-spots of the Fermi surface are ab- 
sent. The eigenvalue corresponding to the AFM instabil- 
ity is decreasing with decreasing temperature at low T, 
so that this instability is again not favored in the ground 
state. The corresponding eigenfunction has its maximum 
at the diagonal of the Brillouin zone. Both eigenvalues 
corresponding to the singlet- and triplet superconduct- 
ing pairing are smaller than in the case of /i = 0.80£, 
although the singlet superconductivity remains the lead- 
ing instability in the T — > limit. The shapes of the 
wavefunctions for singlet and triplet pairing are similar 
to those at /i = 0.80i. The maximum of the wavefunction 
in the singlet superconducting channel is further shifted 
towards the (n, 0) point of the Brillouin zone (i.e. in the 
direction opposite to the hot-spots), while maxima of the 
eigenfunction in the triplet pairing channel arc shifted to- 
wards the diagonal. 

Now we discuss the physical origin of the shape of the 
gap functions found in the numerical analysis. To this 
end, we plot the effective pairing interaction Tj p (k, k') 
in the singlet- and triplet channel as a function of the 
momenta k, k' on the Fermi surface (Fig. 4). One can 
see that at low electron doping (/i = 0.8t) the max- 
ima of the attractive interaction in both, the singlet and 
the triplet channels are located near the hot-spots (Fig. 
4a, b), which leads to the maximum of the pairing gap 
near these points. The difference in 1 patch between the 
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FIG. 4. Contour plots of the pairing interaction rj p (k, k') 
in singlet (a,c) and triplet (b,d) channels at t' = 0.3i, 
U — 3.5t, (i — 0.8t (a,b) and /i = lAt (c,d) as a function 
of the position of the momenta k, k' at the Fermi surface (la- 
beled by the number of patch). Black color corresponds to 
the attractive (negative), white - to the repulsive (positive) 
interaction. White contours show regions with the strongest 
attraction, dashed lines mark the points with k x ,y = — k^ >x . 

position of the hot-spots and the maxima of the gaps 
can be explained in this case by small incommensurabil- 
ity of spin fluctuations and the contribution of the other 
channels of electronic scattering. With increasing doping 
the maximum of the interaction in the singlet channel 
spreads in a broader momentum range (Fig. 4c), which, 
however, leads to almost the same position of the maxima 
of the gap, as for small electron doping. The maximum 
of the attraction in the triplet channel at the points with 
kx.t, = — x shifts towards the diagonal of the Brillouin 
zone, where it is compensated, however, by the strong 
repulsion, which arises at k = k' (Fig. 4d). As a "com- 
promise" , the maxima of the gap are located in the parts 
of the momentum space where neither repulsive nor at- 
tractive interaction is strong, i.e. again remain almost 
unchanged with respect to small doping. 

Therefore, while at small doping the shape of the gap 
functions are determined by the spin fluctuations with 
the wavevector close to Q, at larger dopings the spin 
fluctuations with broader range of momenta start to play 
an important role. Note that the maximum of the singlet 
gap is located away from the Fermi surface points closest 
to (it, 0) and (0, 7r) only for the fillings n > 1; on the hole 
doped side (i.e. at fillings n < 1) the nonmonotonicity of 
the singlet gap function quickly disappears [18]. 

In conclusion, we have investigated the symmetry of 



the leading instabilities and the shape of the correspond- 
ing wavefunctions of the 2D electron doped superconduc- 
tors within the 2D t-t' Hubbard model using as a novel 
tool the combination of the Bethe-Salpeter equation and 
the fRG approach. At U — 3.5t we have found the <i-wave 
pairing instability to be the strongest instability for dop- 
ings S > 0.05 (one can expect that the range of dopings 
where the AFM instability is the leading one, increases 
further with increasing interaction strength). The angu- 
lar dependence of the singlet pairing gap function is found 
to be non-monotonic and its shape is strikingly similar to 
the recent experimental data. The position of the maxi- 
mum of the gap is close to the position of the hot spots at 
low doping level and deviates from it at higher dopings, 
where the antifcrromagnetic fluctuations with wavevec- 
tors with broad momenta range become essential. The 
triplet instability has subleading eigenvalues and f-wave 
like form of the wavefunction. The maxima of the wave- 
function are close to the hot-spots at low electron dopings 
and are close to the diagonal for higher dopings. 
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